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PR3FACE 


This  Memorandum,  the  theoretical  portion  of  which  was  completed  in 
1960,  is  Intended  to  serve  as  a  basic  aid  in  studies  of  orbital  mechanics. 
By  presenting  a  relatively  simple  and  straightforward  method  of  optimum- 
and  nonoptimum-orbit  computation,  it  sho\xld  be  of  use  to  the  engineer  con¬ 
cerned  with  preliminary  mission  analysis.  In  addition,  the  method  of  solu¬ 
tion  described  here  has  pedagogical  merit  because  it  enables  beginning  stu¬ 
dents  of  orbital  mechanics  to  solve  assignments  involving  three-dimensional 
transfer  orbits,  without  the  need  for  machine  computations. 


-V- 


SUMMARy 

The  determination  of  possible  interplanetary  transfer  trajectories  for 
any  specified  date  of  departure,  in  a  realistic  planetai:/  model,  is  a  complex 
problem  'Vfliiah  requires  the  use  of  hl^-speed  ces^ters.  Reference  2  presents 
sooB  of  the  solutions  obtained  for  the  case  of  Mfixtian  transfers. 

For  purposes  of  preliminary  design  and  mission  planning  it  vould  be 
useful  to  be  able  to  compute  some  realistic  transfer  curves  without  extensive 
machine  runs. 

This  Memorandum  presents  a  method  of  solution  which  enables  one  to  de¬ 
termine  arbitrary  transfer  orbits  quickly  and  accurately  using  only  a  desk 
coDguter. 

Coplaoar,  as  veil  as  three-dimensional,  surface-to-surface  trips  were 
investigated,  and  account  was  taken  of  the  eccentric  shape  of  the  planetary 
orbits. 

A  numiber  of  transfer  orbits  from  Earth  to  Mars  vere  con^ted  for  two 
arbitrarily  selected  dates  of  departure.  Ihe  numerical  results  are  sum¬ 
marized  in  a  series  of  curves  which  display  the  characteristic  velocl'ty  ex¬ 
penditures  for  trips  of  various  dixratlons,  as  well,  as  the  orientation  of  the 
departure  velocity  vector. 

In  the  case  of  three-dimensional  transfers,  two  separate  regions  of 
minimal  characteristic  velocity  are  found.  The  two  regions  are  separated 
by  a  narrow  belt  corresponding  to  transfer  angles  in  the  vicinity  of  l80 
deg.  Within  this  zone  of  separation,  which  coiTesponds  to  departures  at 
right  angles  to  the  ecliptic,  transfer  velocities  rise  sharply. 


-vli- 


COMTEHTS 

PKEIFACS .  iii 

SUMMARY . V 

LIST  OF  SYMBOLS . ix 

Section 

I.  INTRODUCTIOT .  1 

II.  GEOMETRIC  PROPERTIES  .  4 

III.  COELANAR  TRAMSFERS  BETWEEN  CIRCULAR  ORBITS .  13 

IV.  THREE-DIMENSIONAL  TRANSFERS  WITH  ECCENTRICITY .  24 

Appendix 

A.  NUMERICAL  EXAMPLE . 35 

B.  CHARACTERISTIC  TRANSFER  VELOCITY  IN  OSIEEE  DIMENSIONS  . .  41 

REFERENCES . 45 


-ix- 


SYMBOLS* 


i 

L,  t 

M 

M,, 

n 

P 

Q 

r 

*'a 

S,  8 


semlnajor  axis  of  elliptlo  orbit 

specif lo  values  of  the  senisBjor  axis,  as  defined  in  Eq.  (3} 

chord  Joining  departure  point  P  to  destination  point  Q 

eecentilcity  of  elliptic  orbit 

active  focus,  or  function  defined  by  Eq«  (lO) 

vacant  focus  located  on  branch  of  Fig,  2 
vacant  focus  located  on  L_  branch  of  Fig.  2 

function  defined  by  E.,  (9) 

included  angle  between  planetary  orbital  planes 

senllatus  rectum 

mean  anomaly  of  destination  planet 

values  of  mean  anomaly  at  departure  and  arrival,  xespeetively 

ascending  and  descending  nodal  axes 

mean  angular  velocity  of  Mars  (NtP. ) 

number  of  complete  revolutions  of  transfer  vehicle 

designation  of  departure  point 

designation  of  arrival  point 

mean  radii  of  planets  Earth  and  Mars,  respectively 
radial  distance  from  Sun 

mean  radii  of  Initial  and  terminal  planetary  orbits 
parameter  defined  after  Eqs.  (l) 


- IF - 

In  the  following,  whenever  a  letter  symbol  appears  as  both  capital 
and  lower-case  (e.g.,  A,  a),  the  capital  refers  to  coodiaenslonallzed 
quantities.  In  all  other  Instances,  nondimensional  parameters  will  be 
Indicated  by  (N.D. ), 
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t  «  tine 
T  ■  tine  (N.D,  ) 

Tjj  =  time  talsen  Tjjt  Itos  to  co'-'er  prescribed  arc  (N.D. ) 

V,  V  «  velocity 
^GH  =  Characteristic  velocity 

■  characteristic  velocity  eoa^ponents  at  deperture  and  arrival 

2 

a,  p  ■  peraaeters  defined  after  Eqs.  (1) 

y  m  angle  subtended  by  ecllx>tlc  con^nent  of  departure 
velocity  vector  vith  Earth  *8  orbit 

8  «  angle  subtended  by  departure  velocity  with  ecliptic  plane 
e  «  angle  included  between  transfer  and  destination  planes 

9  a  transfer  angle 

«  initial  longitude  difference  of  planets  (two-dlnenelonal) 

X  «  constant  lagrange  multiplier 
n  »  solar  gravitational  constant 
V  =  true  anansly 

p  ■  radial  distance  from  8un  (N.D, ) 

Pg  »  mean  radius  of  Martian  orbit 
Pj^  •  radial  distance  corresponding  to  an  anomaly  M 
0  a  final  position  of  Mars  from  ascending  nodal  axis 
*»  angle  traversed  by  Mars,  as  given  in  Eq.  (4) 

■  initial  position  of  tfars 

(■  «  location  of  Mars'  perihelion,  as  shown  in  Pig.  7 
0)  a  ;positlon  of  Earth  at  departure 
«  mean  angular  velocity  of  Mars 


SUBSCRIPTS 
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a  »  aphelion 
CH  ■  Characteristic 
E,  M  ■  Earth  and  Mars^  respectively 
1  «  initial 

FQ  ■  novement  on  transfer  ellipse  from  point  F  to  point  Q 
p  •  perihelion 

Tf  0  »  radial  and  transverse  directions  in  Martian  pImm 

^  # 

X  ■  refers  to  ary  one  of  the  subscripts  F,  P  ,  or  FP  , 
depending  on  circumstances 
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1,  INTRODUCTION 

During  the  planning  stage  of  any  interplanetary-transfer  miasion,  there 
exists,  for  purposes  of  prelimlnaxy  design,  a  very  distinct  need  for  quich 
and  reliable  information  concerning  the  merits  of  the  various  transfer  orbits 
'which  can  be  selected  to  accocpUsh  the  given  mission. 

Two  of  the  main  items  of  interest  during  this  stage  of  the  investigation 
are  the  time  of  fli^t  and  -the  so-called  characteristic  velocity  of  transfer. 

Tine,  ob'viousiy,  is  of  great  importance  if  one  considers  manned  space 
miasiona,  since  it  affects  directiy  the  crew  requirements  for  the  duration 
of  the  trip. 

The  characteristic  velocity,  on  the  other  hand,  affects  the  propellant 
requirements,  and  thus  the  size  of  the  final  vehicle.  The  influence  on 
vehicle  size  decreases  markedly  as  one  progresses  to  the  use  of  fuels  having 
a  vary  high  specific  impulse. 

In  the  Initial  stages  of  interplanetary  flight  analysis,  mission-cost 
and  flliidit-tlme  estimates  between  coplanar,  concentric  circular  plane'tazy 
orbits  were  obtained  by  the  application  of  the  concept  of  cotangentlal  transfer 
orbits,  also  known  as  Hohmann  ellipses.  Althot;^  these  transfer  orbits  led 
to  low  ■values  of  the  characteristic  transfer  velocity,  they  required  rela¬ 
tively  long  ti*an3fer  times.  In  addition,  the  use  of  thes«i  orbits  presupposes 
a  propitious  location  of  -the  t-wo  planets  at  departure,  a  condition  seldom 
encountered  in  practice. 

Cotangentlal  ellipses  cease  to  exist  in  a  three-dimensional  model  of  -the 
planetary  sys'tam,  irrespecti-ve  of  -whether  the  planetary  orbits  are  still  con¬ 
sidered  to  be  circular  or  elliptic.  Even  so,  this  concept  can  still  find  an 


occasional  application  for  those  Instances  \Aien  the  target  planet  happens 
to  cross  the  nodal  axis  at  the  time  of  Intercept,  and  the  departure  point 
is  sltmted  l8o°  away. 

With  the  exception  of  these  special  cases,  transfer  orbits  will  lie  In 
planes  of  arhltraiy,  though  determinable,  oxdentatlons  In  space,  depending 
on  the  departure  date  selected  and  the  time  available  to  transit.  Ibe  simple 
bach-of-the-envelope  type  of  trajectory  comf/utations  for  the  two-dimensional 
circular  planetary  model  that  were  Introduced  by  the  early  Investigations  in 
the  field  were  replaced  by  cosq^cated  coding  routines  for  hl^-speed  digital 
coDQuters,  Into  idilch  the  exact  orbital  positions  of  the  planets  could  be 
fed. 

Ihis  Memorandum  presents  a  method  vdilch  allows  one  to  determine  and 
cooipute  with  relative  ease  any  type  of  two-impulse  transfer  trajectory,  in 
two  and  three  dimensions,  for  a  given  date  of  departure  if  the  time  avail¬ 
able  for  the  flight  Is  specified  In  advance.  Since  only  elliptic  transfer 
orbits  have  been  considered  here,  this  iioplles  that  among  all  the  possible 
ellipses  which  pass  throu^  the  given  Initial  point,  and  specified  final 
point,  only  that  ellipse  is  Iscilated  which  enables  the  vehicle  to  spend  the 
pre;-clected  nusher  of  days  in  transit.  Moreover,  this  calculation  can  be 
carried  out  In  a  very  short  period  of  time  (2  to  3  hours)  on  a  regular*  desk 
computer  without  the  use  of  high-speed  digital  machines.  In  addition,  it 
Is  shown  how  the  formallsn  of  some  of  the  theorems  of  the  calculus  of  ordi¬ 
nary  maxima  and  mlniina  facilitates  the  location  of  those  trajectories  which, 
for  a  given  departure  date,  require  the  least  expenditure  of  characteristic 
velocity  of  transfer. 

In  Ref.  1,  Battin  introduced  a  very  convenient  geanetrlc  formulation  in 
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the  two-diiaensional  portion  of  his  irork,  in  -iAich  he  treated  transfers  be- 
tveen  a3*itrary  points  on  circular  planetary  orbits.  She  necessary  initial 
locations  of  the  planets  were  then  computed  as  soon  as  an  optimal  trajectory 
'uas  found.  Ihe  same  geometric  parameters  will  be  retained  here. 

Starting  out  with  a  two-dimensional  problem,  we  introduce  a  subsidiary 
intercept  constraint  condition  on  the  transfer  orbit  vhich  enables  one  to 
calculate  transfers  directly  for  only  those  planetary  configurations  at  de¬ 
parture  ■vliich  are  of  current  interest.  The  formalism  of  the  optimality  con¬ 
dition  for  the  present  case,  analogous  to  the  one  developed  by  Battin,  is 
some^Aat  more  complicated  because  of  the  need  to  Introduce  Lagrange  multi¬ 
pliers. 

13re  same  basic  approach  is  then  extended  to  solve  the  case  vhere  the 
destination  orbit  is  elliptic  in  shape,  and  it  is  eventually  applied  to  the 
more  realistic  model  in  lAiich  the  orbital  planes  of  the  planets  are  naituoUy 
Inclined  to  each  other. 

Qie  difficulties  introduced  by  tlie  Inclusion  of  orbital  eccentrlcily 
have  been  'overccxrae  by  retaining  in  the  series  exsansions  of  the  planetary 
orbital  elements  only  first-order  terms  in  eccentricity.  Ihe  very  small 
valiies  of  this  quantity  for  the  orbits  of  Earth  and  Mars  Justity  this  approxl 
matlon. 

Ihe  numerical  results  of  the  present  investigation  were  compared  with 
those  of  Ref.  ?  and,  to  the  order  of  accuracy  of  the  data  presented  there, 
were  found  to  be  in  good  agreement  with  them. 
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geametrlcal  relationships  'v&lc^  will  be  of  use  In  the  subsequent 
develppnent  vUl  be  stated  here  vlthout  proof.  Ihose  interested  in  the  deri¬ 
vations  are  referred  to  Ref.  1.  Distances,  velocities,  and  times  are  non- 
dimenalonallzed  vith  respect  to  the  astronomical  unit;  Eartii's  mean  speed 
around  the  sun,  ^  }  and  its  mean  angular  velocity  /  ^  »  respectively. 

V  ^1  V  r3 

For  that  part  of  the  calculations  in  Tjhich  the  orbits  of  both  Earth  and 
itos  vere  assumed  to  be  circular,  a  transfer  from  a  point  P  on  the  orbit  of 
the  first  planet  to  a  point  Q  on  the  orbit  of  the  second  planeP  vould  take 
place  between  the  radial  distances  r^^  and  Tg  of  the  transfer  ellipse.  Only 
direct  transfer  orbits  vere  considered  here  because  of  their  lover  cost. 

For  a  prescribed  angle  0  between  the  above  two  radii,  the  time  of  transit 
will  depend  on  both  the  magnitude  and  the  orientation  of  the  semlmajor  axis 
of  the  transfer  ellipse.  (Hie  time  can  be  calculated  from  Lambert's  theorem, 
a  derivation  of  lAiich  can  be  found  in  Ref.  3.  Denoting  by  F  the  active  focus 
of  the  transfer  ellipse  at  which  the  Sun  is  located,  and  by  F  the  vacant 
focus,  we  can  dlstlngolsh  among  four  possible  modes  of  transfer  from  P  to  Q, 
as  indicated  in  Fig.  1. 

Which  one  of  the  following  relations  for  the  time  of  travel  will  be  used 
at  any  one  instance  will  depend  on  i^iether  the  contour  from  P  to  Q  and  back 
to  P  via  the  chord  QP,  traversed  in  a  counterclockwise  direction,  does  or 
does  not  enclose  one  or  both  of  the  focal  points  F  and  F  , 

Hie  sense  of  the  motion  will  be  taken  as  positive  in  the  counterclock¬ 


wise  direction. 


Hie  various  expressions  for  the  dimensionless  time  of  transit  can  then 
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be  written  as  follows: 

T  -  a3/2 

^PQ  * 


(a  -  sin  Qt)  -  (p  -  sin  p) 

Tpp*Q  »  2nA3/2  _  ^3/2  .  gin  j^)  +  (3  _  g^n  p) 


■^EFQ 


,3/2 


(a  -  sin  a  +  (p  -  sin  p 


"  2nA3/2  . 


where  the  angles  a  and  p  are  given  by  the  relations 


I  aa 


§ .  s.s^ 


2a  ■  r^  +  rg  +  c 


T  »  t 

^1 


A  & 


(la) 

(lb) 
(lo) 
(Id) 


p  <  a  <  TT 

Expression  (la)  is  used  %^en  the  contour  does  not  enclose  either  one  of  the 
fooi|  as  in  (a)  of  Fig.  1.  VJhen  FQF  encloses  the  vacant  focus  F  ,  as  in  (b) 
of  Fig.  1,  one  uses  Expression  (lb).  IS,  on  the  other  hand>  the  aoti^re  focus 
is  enclosadj,  as  in  (0)  of  Fig.  1,  one  uses  Expression  (lo).  VOien  both  F  and 
F  are  situated  within  the  contour  F<2F,  as  in  (d)  of  Fig.  1,  one  oust  apply 
Expression  (Id). 

Another  useful  paratneter  is  the  semllatus  rectua  L  of  the  transfer  orbit. 
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■where  ■  perifocal  radius,  and  the  an^es  a  and  g  are  functions  of 

Tl  +  r. 

■  •  ■  ^  and  C. 

For  thoae  oases  idiere  the  tiiae  of  transit  did  not  have  to  satisfV  a 
constraint  condition,  the  selection  of  the  positive  sign  in  the  e3q)re88ion 
for  the  aemilatus  rectxan  L  led  to  lower  values  of  the  cJmracteristics  trans¬ 
fer  velocity.  Hhis  can  be  shown  by  plotting  vs  A,  and  it  is  found  to 
hold  true  for  all  values  of  the  transfer  angle  0.  Now  that  the  transfer 
time  has  a  constraint  in^sed  on  it,  it  can  be  shown  that  not  8.11  the  dates 
considered  for  departure  allow  one  to  follow  an  elliptic  path  with  L^.  To 
ensure  the  existence  of  elliptic  transfer  orbits  for  all  the  possible  initial 
positions  of  the  planets,  and  to  minimize  the  number  of  cos^ete  orbital  revo¬ 
lutions  before  arrival,  both  valiies  of  the  quantily  L  are  retained  in  the 
analysis. 

One  further  point  worth  mentioning  concerns  the  orientation  of  the  trans¬ 
fer  ellipse.  For  given  values  of  r^^,  r^,  and  included  angle  0,  this  orien¬ 
tation  depends  on  the  magnitude  of  the  semi  major  axis,  as  well  as  on  the  po- 

♦ 

sitlon  of  the  •vacant  focus  F  •wl'th  respect  to  -the  chord  C  Joining  P  and  Q. 

■Die  choice  of  the  sign  in  the  eigireosion  for  the  seiiil.latus  rectum  serves  to 
determine  on  which  side  of  the  chord  C  the  focus  P  will  be  looa-ted.  Figure  2, 
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>Ailch  1b  similar  to  Figs.  1  -  4  of  Ref.  1,  indicates  the  locus  of  possible 
* 

F  positions,  vlth  the  heavy  bxanhh  of  the  curve  denoting  the  region  of 

nils  curve  turns  out  to  be  a  hyperbola.  If  a  particular  position  of  P  is 

chosen,  the  orientation  of  the  major  axis  of  the  transfer  ellipse  is  then 

* 

found  by  Joining  the  points  F  and  F  .  Depending  on  the  orientation  of  this 
axis,  a  vehicle  leaving  point  P  along  a  suitable  transfer  orbit  vould  reach 
point  Q  either  before  or  after  passing  an  apse  of  the  transfer  ellipse. 

Referring  to  Fig.  2,  a  fev  critical  points  located  on  the  curve  de> 
temined  by  the  loci  of  F  can  be  noticed.  !]he  point  F  is  the  vacant  focus 

ID 

of  the  ellipse  with  the  minimum  semlmajor  axis;  F^  is  the  vacant  focus  of 
the  transfer  ellipse  with  aphelion  located  at  point  Q.  For  values  of  ^  <  6  < 
a  third  point  of  interest  is  F^^  vhleli  defines  'die  ellipse  with  perihelion 
at  P.  With  each  of  these  vacant  foci  there  is  associated  a  magnitude  of  the 
semlnajor  axis  of  the  transfer  ellipse.  Denoting  these  respective  parameters 
by  A^,  and  A^,  where 


T,  +  r-  +  c 
- 

_ _ 


1  + 


rg  -  cos  9 


r^  -  Tg  cos  9 

2r^  -  rg(l  +  COB  0 ) 


.31! 


(3) 


and 


m 
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ve  can  cletermine  the  orientation  of  the  transfer  e^Jdpse  from  P  to  Q  sisiply 
from  a  knowledge  of  the  magnittide  of  its  senimajor  axis  A. 

The  dependence  of  travel  time  t  on  semimajor  axis  A  of  the  transfer 
ellipse  is  indicated  in  Fig,  3  for  two  arbitrajry  values  of  the  transfer 
angle  6. 

Turning  our  attention  for  the  moment  back  to  Pig,  a,  we  observe  that 

* 

for  a  given  x^alue  of  8,  the  junction  of  the  two  branches  of  the  loci  of  F 

occurs  at  a  value  of  A  »  A_.  ,  The  aanw  behavior  is  also  exhibited  in  the 

jmn 

X  vs  A  graph  of  Fig,  3* 

From  Fig,  2,  -vdiich  happens  to  be  drawn  for  9  >  ■we  note  that  as 

* 

A  -)<»,  the  direction  of  the  major  axis,  or  the  ll.ne  joining  P  and  F  on  the 

* 

branch,  becomes  parallel  to  the  asyniptote  of  the  hyperbola  of  F  ,  This 
places  the  perifocal  apse  on  the  segment  of  arc  Joining  P  and  Q.  Since  this 
region  implies  hi^  orbital  velocities,  it  is  reasonable  to  assume  that 
transfer  times  along  this  arc  are  going  to  be  of  short  duration.  This  ex¬ 
pectation  is  borne  out  by  the  branch  of  the  9  «  60  deg  curve  of  Fig.  3. 
(The  fact  that  here  9  <  ^  will  not  alter  the  conclusions  reached  before, 
except  for  the  fact  that  now  the  perifocal  apse  will  be  located  slightly 
oliead  of  the  point  P,  rather  than  between  P  and  Q.  Point  P  will  still  have 
hi(^  orbital  velocities  associated  with  it. )  As  the  point  F  slides  along 
the  branch  of  the  hyperbola  in  the  direction  of  decreasing  A,  the  transfer 
ellipses  rotate  in  a  clockwise  direction,  thereby  moving  the  perifocal  apse 
away  from  point  P.  The  attendant  lower  orbital  velocities  in  the  region 
between  P  and  Q  vrill  result  In  on  increase  in  travel  time  t,  as  exemplified 
by  the  lower  branch  of  the  0  ■  60  deg  cxirve.  The  largest  available  value 
of  -v  along  this  branch  will  octnur  when  a  value  of  .A  ■  A 


is  reached.  If 


11. 


Fig.  3  —  Variation  of  transfer  time  with  semimajor  axis 
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longer  travel  tinea  are  desired,  the  vacant  focus  P  -will  have  to  continue 
to  slide  out  into  the  P  region  (L  »  L_),  thereby  causing  the  vehicle  to 
move  past  a  more  remote  apofocal  point,  and.  conseoiuently  into  a  region  of 
ever-diminiahing  orbital  velocities.  !Qie  lover  orbital  velocities  coupled, 
vith  the  longer  distances  that  have  to  be  traversed  as  k  keeps  on  increasing 
will  lead  to  longer  transfer  tlaea,  as  shown  in  the  upper  branch  (L_)  of  the 
0  =  60  deg  curve  of  Fig.  3. 

For  values  of  0  >  tt,  the  and  L_  branches  of  the  hyperbola  become 
reversed,  thus  explaining  the  6  =  240  deg  curve  of  Fig.  3. 

An  examination  of  the  equations  presented  in  this  section  discloses 
the  fact  that  for  transfers  between  circular  orbits  where  and  r^  are 
constants,  the  time  t  as  well  as  the  semllatus  rectum  L  depends  only  on  the 
semimajor  axis  A  and  the  chord  C  between  points  P  and  Q. 
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ni.  COPIi/mAR  TRAMSFBRS  BETWEEN  CIRCULAB  ORBITS 

Before  looking  at  the  more  general  and  realistic  case  of  transfers  In 
three  dimensions  betveen  elliptic  orbits,  It  Is  advantageous  to  start  out  tgr 
considering  at  first  the  simpler  case  of  transfer  from  the  surface  of  the 
Barth  to  the  surface  of  Mars,  assuming  that  both  planets  describe  circular 
coplanar  orbits  about  the  Sun.  With  some  obvious  minor  modifications,  the 
same  analysis  would  apply  also  to  transfers  to  inner  planets,  such  as  from 
Barth  to  Venus,  or  to  tfc.e  return  trip  from  Mars  to  Earth. 

It  vUl  be  assumed  that  the  mission  Is  accon^lished  by  the  application 
of  two  Inpilslve  thrusts  at  both  termlnala  of  the  transfer  curve.  Qie  first 
thrust  enables  the  v^lcle  to  escape  the  Earth’s  gravitational  field  and 
leaves  it  with  a  hyperbolic  excess  velocily  of  a  magnitude  and  direction  such 
as  to  place  it  Into  heliocentric  elliptic  intercept  orbit.  Ihe  second  Im¬ 
pulse  Is  emplciyed  both  to  cancel  the  relative  velocity  of  the  vehicle  vrlth 
respect  to  Mars,  on  Its  arrival  there  along  the  cdiosen  heliocentric  collision 
course,  and  to  overcame  the  gravitational  attraction  of  that  planet.  Oils 
will  pentlt  the  vehicle  to  make  a  soft  liuidlng  on  Mars.  Ihe  assumed  scheme 
of  transfer  does  not  make  euiy  allowances  for  the  various  losses  ( gravitati<onal, 
atmospheric,  navigational,  etc. )  that  arise  in  an  actual  space  mission  nor 
for  possible  Egains  from  atmospheric  braking  on  re-entxy.  Nonetheless,  the 
present  method  will  be  well  suited  to  assess  the  rou^  over-all  costs  of  sudi 
a  transfer  and  to  give  an  idea  of  the  correlation  between  launch  date  and 
total  velocity  requli'«ment.  hold  the  needed  characteristic  velocities 
down  to  a  reasonable  Itivel,  transfers  along  elliptic  orbits  only  will  be  In¬ 
vestigated.  dhe  modifications  needed  to  include  faster  transfers  will  affect 
pzimarily  the  expressions  for  time  of  transfer  and  semllatus  rectum. 
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Ihe  complete  siirface-to-surface  transfer  trajectory  tos  assumed  to  be 
coursed  of  three  discrete  segments,  each  segment  being  generated  as  a  re¬ 
sult  of  a  tvo-body  interaction. 

The  junction  of  the  heliocentric  portion  of  the  orbit  with  the  planet- 
ocentrlc  segments  at  departure  aM  arrival  was  assumed  to  occur  along  the 
asymptotes  to  the  planetooentric  hyperbolae.  At  each  end,  the  planetary 
influence  on  the  terminal  portion  of  the  trajectory  was  assumed  to  be  con¬ 
centrated  in  an  imposition  of  the  appropriate  escape  velocity  requirement, 
■vAilch  was  added  vectorialljr  to  the  planetooentric  hyperbolic  excess  velocity. 

'The  major  benefit  vAilch  the  present  approach  brings  to  the  three-di¬ 
mensional  case  is  the  ability  to  solve  for  any  one  transfer  trajectory,  given 
the  positions  of  the  two  planets  at  departure. 

Consider  the  physical  picture  -iftiere  the  two  planets  occupy  specified 
positions  along  their  orbits  on  the  specific  date  chosen  for  launch  from  the 
Earth.  In  the  two-dimensional  restricted  case,  only  the  angle  between  the 
orbltel  radii  of  the  two  planets  suffices  to  specify  their  relative  position 
at  departure,  their  actual  positions  in  their  respective  orbits  being  im¬ 
material. 

After  the  initial  impulse  is  imparted  to  the  vehicle,  it  will  find  it¬ 
self  with  a  certain  initial  velocity  vector,  in  a  heliocentric  frsaie  of  refer¬ 
ence.  IMs  velocity  is  already  assumed  to  include  the  contribution  from  -the 
circumferential  velocity  of  the  Earth  in  its  path  around  the  Sun. 

The  two  components  of  this  departure  velocity  vector,  in  the  radial  and 
transverse  directions,  suffice  to  specify  completely  the  geometrical  para¬ 
meters  of  the  transfer  ellipse.  Starting  out  with  the  above  conditions  at 
point  P,  and  following  the  selected  orbit,  the  vehicle  will  intersect  the 
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orMt  of  Mars  at  sotas  future  time,  at  Boes  point  Q.  Since  this  crossing 

aust  occur  at  the  precise  ament  vhen  Mars  happens  to  pass  the  saaie  spot,  it 

hecooes  neeesaaiy  to  restrict  the  domain  of  posslhlc  departure  conditions. 

A  sufficient  restriction  on  the  departure  conditions  results  if  one  demands 

that  the  tiaie  of  travel  from  P  to  along  the  transfer  orbit,  be  the  sane 

as  the  tins  required  by  Mars  to  reach  the  saae  terminal  point  Q*  If  one  as- 

soclates  vlth  each  tine  of  transfer,  defined  by  the  appropriate  eq^aallty  in 

Bqs.  (l),  a  transfer  angle  0  included  betveen  the  radii  r,  »  r„  and.  r_  «  r„ 

X  X*  2  H 

and  remembers  that  in  the  present  case  Mars  moves  with  a  constant  angular 
veloci1;y,  one  can  then  vrlto 


vhere 


^  m  Q  ^  0  -  arm 
o 

Vj. 

a  ~  t  •  angle  traversed  by  >farB  during  time  t 
2 

«  lead  angle  of  ^fars  at  departure 

Vjj  ■  dimensional  circuirferential  velocity  of  Mars 

n  •»  0,  1,  2,  3  "  nuuiber  of  revolutions  around  the 
Sun  completed  by  the  vehicle 


(>^) 


For  the  selection  of  optimal  transfer  curves,  the  proced'ure  of  optimization 
then  consists  of  hunting  among  the  various  acceptable  points  Q  for  tJiat 
particular  point  vhlch  permltB  the  mission  to  be  acconqtUshed  vlth  a  ndnlimnu 
e}(penditure  of  total  characteristic  velocity.  IMs  hunting  csui  be  performed 
airtomatlcally  by  means  of  a  hl^-speed  computer,  or  manually  on  a  desk  calcu¬ 


lator, 
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Another  possitle  vay  of  viewing  the  problem  is  the  one  in  'vAiich  we 
seardi  for  that  particular  transfer  ellipse  v*idi  requires  the  vehicle  to 
spend  a  period  of  time  in  orbit  of  such  length  that  both  the  requirement  of 
minimal  <aiaraoteristic  velocity  and  the  restraint  imposed  by  Eq.  (l^)  are 
satisfied  simultaneously. 

With  the  above  physical  picture  in  mind,  it  is  easy  to  appreciate  the 
mathematical  fonmilatlon  which  leads  to  the  required  solution. 

The  quantity  to  be  minimized  in  the  present  case  is  the  total  cliarac- 
teristio  velocity  of  the  one-way  landing  mission,  given  in  nondimensional 
form  by 


V  {3  -  i  -  sVi  *  (5) 


^1 

a  -j^  ^  a  square  of  dimensionless  escape  velocity  at  Earth 

2uj^  ^1 

K,  =  — —  ^  »  square  of  dimensionless  escape  velocity  at  Mars 

4 


The  first  term  of  Eq.  (5)  was  the  one  derived  in  Ref.  1  for  a  nonattraotlng 
Earth  (Kg  »  0). 

The  chord  joining  points  P  and  Q  can  be  expressed  in  terras  of  the  angle 
0  by  means  of 


=»  1 


2p^  cos  0 


(6) 


where  the  angle  0  depends  on  the  time  of  travel  of  Mars,  through  the 
relation 


O-2nnB0  +N,,  T 
o  K  o 


(7) 
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In  the  above  eqp»tion,  dlmenBionless  mean  angular  ve- 

locily  of  Mars,  and  is  the  corresponding  dimensional  angular  -velocity. 

Ihe  reference  time,  la  58,12  days. 

Ihe  selection  of  a  departure  da-be  fixes  Immedia-tely  -the  -value  of  -ilte 
angle  9^.  In  order  -bo  mahe  contact  vith  Mars  on  arri-val,  one  has  to  reqtiire 
that  -the  time  of  travel  in  the  -transfer  orbit,  t,  given  by  -the  appropria-te 
e^^nresslan  in  Eg,B.  (1),  should  equal  -the  time  -taken  by  Mara  -to  reach  the 
same  spot,  for  e-vexy  value  of  -the  angle  6^. 

Qius 

T  -  Tjj  (8) 

is  a  required  condition  on  -the  transfer  orbit.  Noting  -that  by  maana  of 
Eqs,  (6)  and  {7}  ~  f^rom  Sqs,  (l)  the  time 

T  ■  t(A,C),  one  can  -Hrite  CJondition  (8)  in  -the  form 

s(A,C,0^)  -  T  -  Tjj  »  0  (9) 

Our  optimization  problem  nov  -takes  -the  faULowlng  form. 

It  is  desired  -to  mlntmlza  -the  function  ■  Vjjjj(A,C)  given  in  Eq.  (5), 
subject  to  -the  subsldiazy  condition  of  Kq.  (9}.  In  accordance  vith  a-tandard 
pirooedures  of  the  Siculus  of  max! ns  and  minima,  this  entails  the  determi¬ 
nation  of  an  extremal  value  of  a  function  E,  given  by 

F  -  V(jH  +  Ig  (10) 

vbere  X  Is  a  constant  Lagrange  multiplier.  Ihe  extremals  are  obtained  upon 
setting 


8p8f8F- 
3a  "'3c  ° 


(11) 


-18- 


The  relation 


&F  ^'^CH  ,  .  ^ 

5c  “  "5c"  5c 


can  1)6  used  to  sol'V'e  for  vihlch  can  then  be  substituted  in  the  equation 

Sf 

^  B  0,  l-lhen  this  Is  done^  the  condition  for  an  extremal  in  becames 


(12) 


in  addition  to  Qq.  (9). 

A  suitable  method  for  computing  optimal  trajectories  -would  consist  of 
-the  following  s-teps: 

a.  For  a  particular  date  of  departure  the  angle  is  a  known  quantity, 
a-vailable  from  astronomical  tables.  Ihe  angular  positions  of  Earth  and  Mars, 
reckoned  from  -the  ascending  nodal  line,  are  prescn-ted  in  Fig.  4>  In  -the  case 
of  circular  orbits,  both  ctirves  become  stral^t  lines  of  known  slope.  IHiese 
lines  Intersec-ted  on  January  1,  19^1,  at  which  time  Mars  was  in  opposition 

to  the  Earth. 

b.  A  guess  for  a  reasonable  -value  of  transfer  angle  9  allows  one  to 
coD^JU-te  C  and  by  means  of  Eqa.  (6)  and  (7),  respecti-vely. 

c.  In  -view  of  Eq.  (9),  Eq.  (l)  can  then  be  solved  for  a  -value  of  the 
semlmajor  axis  A.  nils  s-tep  will  In  general  In-vol-ve  an  Iteration  procedure 
for  A,  -which  con-verges  rather  rapidly. 

d.  The  -values  of  A  and  C  found  from  steps  (b)  and  (e)  above  are  then 
placed  Into  Eq.  (12),  and  the  -value  of  T(A,C)  is  coc^red  wl-th  zero, 

e.  the  procedure  is  repeated,  starting  cTut  wi-th  a  sll^-Uy  changed 
value  of  6,  while  keeping  0^  unchanged. 


Fig.4— Position  of  Earth  and  Mars  referred  to  ascending  node 


f,  graiii  obtaiiied  by  plotting  T(A,C)  vs  Q  waB  found  to  be  well 
behaved;  and  as  a  consequence  a  tedinlque  like  the  Nevton-Baphson  method  of 
successive  approximations  could  be  utilized  to  speed  up  the  convergence  of 
the  solution  X(d)  «  o.  (The  dependence  of  T  on  d  arises  after  the  solution 
of  the  time  £!q<  (l)  for  A  has  been  carried  out;  so  that  A  Itself  la  a  function 
of  0.) 

Where  one  is  not  primarily  interested  In  finding  only  optimal  modes  of 
transfer;  but  desires  to  get  a  fuller  picture  of  the  various  trade-off  values 
resulting  from  changes  In  Initial  conditions;  one  can  rreglect  steps  (d)  - 
(f)  and  compute  ^CH  after  step  (o).  In  this  way  it  Is  possible  to  see  how 
the  charaoterlstlc  veloci'^  varies  with  the  transfer  angle  0.  Figure  ^ 
presents  a  plot  of  against  travel  time  for  two  chosen  departure  dates. 
Figure  6  presents  the  two  con^nnsnts  and  of  the  total  ti'ansfer  re- 
(julrenent:  denotes  the  charactsristio  velocity  at  launch;  while  V^j 

refers  to  the  landing  portion. 

The  orientation  of  the  transfer  ellipse  at  the  point  of  departure  can 
be  deteradnad  with  the  aid  of  the  angle  7^  ^Ich  the  Initial  velocity  vector 
DELkes  with  the  orbit  of  the  Earth,  as  found  from 


tan  ■ 


ihe  calculated  values  of  fix  the  magnitudes  of  the  hyperbolic  excess 
velocities  idiloh  agreed  closely  with  those  of  Fig.  1;  Ref*  2;  within  the 
accuracy  to  ;dilch  the  data  are  presented. 

As  can  be  seen  from  the  preceding  paragraphs;  the  procedure  of  select¬ 


ing  optimal  trajectories  for  the  present  two-dimensional  models  differs  con¬ 
siderably  from  the  one  employed  in  Ref.  1.  Whereas  there,  due  to  the  absence 
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Fig. 5 — Coplanar  characteristic -velocity  variation 


( dimensionless 
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Fig.6  — Coplanar  departure  and  arrival  velocity  components 
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of  a  constraint  on  the  travel  time,  all  the  ellipses  posse a sing  an  A  > 
could  "be  considered  as  suitable  transfer  orbits  for  any  selected  transfer 
angle  9,  ttie  present  specification  of  an  initial  position  of  the  target 
planet  in  general  reduces  this  nunber  to  a  single  suitable  ellipse.  Further- 
more,  in  Bef.  1  it  was  shown  that  there  existed  an  optimal  transfer  ellipse 
for  every  se]«cted  angle  of  transfer.  This  is  no  longer  true  >ftiere  the  seni- 
major  axis  A  is  dependent  on  the  transfer  angle  chosen  through  the  time-con¬ 
straint  condition.  Optimal  modes  of  transfer  can  now  occur  only  for  discrete 
values  of  the  transfer  angle,  or,  vdie,t  amounts  to  the  same  thing,  transfer 
time. 

For  purposes  of  illustrating  the  method  of  solution  in  the  opplanar  case, 
a  numerloal  example  is  worked  out  in  Appendix  A. 
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IV«  TgRBE-DIMEMSlOKAL  ‘TRAySFERS  WITH  BCCBiraRICITy 

to  mm  m  have  eonsidexed  only  circular  coplanar  planstazy  orbits, 
ilhese  restrictions  on  -Oie  orbits  vill  nov  be  lifted,  thereby  extending  con¬ 
siderably  the  scope  of  the  method  developed.  Ibis  vlU  permit  one  to  de¬ 
termine  transfer  trajectories  between  elliptic  noncpplanar  planetary  orbits. 
Hie  elUptle  shape  of  tiie  orbit  of  the  destination  planet  will  affect  the 
manner  of  computing  the  travel  tine  of  that  planet,  while  the  presence 
of  eccentricity  in  the  initial  orbit  will  change  the  form  of  the  e:q>reasion 
for  the  velocity  increnent  to  be  added  at  departure,  V~  .  In  the  majority 
of  eases  the  eccentricities  of  the  planetary  orbits  are  rather  small  quanti¬ 
ties.  Ibis  BBhes  it  possible  to  simplify  natters  appreciabl;/'  by  retaining 
only  linear  terms  in  e  in  the  series  expansions  of  the  planetary  orbital 
elements.  In  particular,  if  one  wants  to  be  consistent  in  the  order  of 
magnitude  of  the  perturbation  terms  retained  for  the  case  of  travel  from 

O 

Earth  to  Mars,  the  fact  that  ="  0  [e^]  allows  one  to  ignore  the  eccen¬ 
tricity  of  the  Barth's  orbit  when  retaining  only  linear  terms  in  the  Martian 
elements.  This  section  will  thus  be  concerned  with  transfers  from  a  point 
situated  on  an  initial  circular  departure  orbit  to  another  point  located  on 
an  inclined  elliptic  destination  orbit. 

Ihe  angle  t  measured  counterclockwise  from  the  ascending  nodal  line  to 
the  perihelion  of  the  destination  orbit  will  be  assumed  to  be  fixed  in  magni¬ 
tude.  For  Mars,  ^  ^  286  deg.  The  geometry  of  the  present  model  is  shown  in 
Figs.  7  and  8. 

In  three-dimensional  maneuvers,  tlie  transfer  plane  will  In  general  sub¬ 
tend  an  angla  5  with  the  ecliptic  and  an  angle  e  with  the  orbital  plane  of 


Fig.7 — Parameters  of  Earth -Mars  model 
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the  destination  planet.  From  sjjherical  trigonometry 

cos  e  «  sin  i  sin  5  cos  ci>  +  cos  i  cos  8  (llj.) 

-vdrere  1  denotes  the  inclination  of  the  destination  plane  vlth  the  ecliptic 
and  Si  denotes  the  location  of  the  Earth  referred  to  the  ascending  nodal 
line. 

Die  transfer  angle  Q  is  now  supplied  by  the  es^ression 

cos  9  =»  cos  <Ji  cos  0  +  sin  cu  sin  0  cos  i  (15) 

in  ■which  0  denotes  -the  angular  position  of  the  destination  planet  at  the 
time  of  intercept,  also  measured  from  the  ascending  node. 

In  ■the  present  case,  due  to  •the  snfill  magni-tude  of  ■the  angle  i,  it  is 
reasonable  to  make  the  approximation  cos  i  1,  after  which  Eq.  (15)  becomes 

cos  0  a  cos  (0  -  to)  (16) 

and  hence 


0  ±  2TTn  ■=  0  - 


(« 


(17) 


Wi-th  -the  mean  anomaly  denoted  by  M,  ■the  dimensionless  radial  distance 
of  Mars  from  the  Sun  becomes 


Pm  "  “  ®M  ” 

“  [2-  -  cos  (t  -  Tp) 


+  0 


(18) 


where 


Pa  Pp 

A,.  «  — r - ^  “  P_  «  radius  of  approximated  circular  orbit 

7u  2  2 
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>■  tins  of  perihelion  passage 
radii  at  aphelion  and  perihelion,  respectively 

w  P 

At  tine  T  a  T^j  vMdh  corresponds  to  departure,  the  true  anomaly  of  Mars, 
vj^.  Is  given  by 

(19) 


tdiere  denotes  the  Initial  mean  anomaly.  If  we  follow  the  convention  that 
positive  angles  are  measured  In  the  counterclockwise  direction,  a  negative 
value  for  would  Imply  that  this  angle  is  aeasxired  clockwise.  In  order  to 
prevent  this  from  happening,  the  lostrlction  (Mod.  Sir)  should  be  understood 
to  naan  that  an  angle  Sir  has  to  be  added  to  the  conputed  value  ■*  i  >Aien- 
ever  this  qiuantlty  becomes  negative.  If  this  is  done,  then  will  always 
be  located  In  the  Interval  0  -  Sv* 

Equation  (19)  Is  basically  Kepler’s  equation  and  can  therefore  be  Inverted 
to  solve  approxinately  for  M^^ 

a  VjL  -  SCjj  sin  VjL  *  (sO) 

Similarly,  at  arrival  t  and  Vg"0-*"(cu-*)+e  -  Sim  (n  »  0,1,2, .. ), 
80  that  for  the  case  n  =  0 


«2 


2e. 


«)  +  e 


(21) 


If  we  understemd  0  to  denote  the  total  emgle  covered  by  the  radius  vector 
from  the  Sun  \dilch  rotates  in  a  counterclockwise  direction,  starting  from 
the  ascending  node  of  Mars  and  stopping  at  the  point  of  Intercept,  after 
having  passed  In  transit  Msrs'  Initial  position  0^,  it  Is  clear  that 
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0  m  0^  +  As  used  here,  A0  represents  the  angle  covered  hy  Mars  during 
the  transfer  time  and  is  a  positive  quantity.  Adhering  to  the  sign  con¬ 
vention  mentioned  in  connection  vith  Eq,  (19),  ve  see  that  vg  >  and  that, 
in  addition,  Vg  can  exceed  2ir  in  magnitude. 

Ihe  time  of  transit  from  the  initial  point  P  to  the  final  point  Q  is 
readily  available  from  the  relation 


exq“  n„ 


(22) 


idiere 


The  letter  X  in  the  subscript  stands  for  any  one  of  the  intenaedlate 
stations  F  or  F  discussed  earlier. 

The  inclination  S  of  the  transfer  plane  is  obtained  from  the  relation 

sin  6  «  sin  i  (23) 

If  ve  refer  to  Eq.  (17),  Eq.  (23)  becomes 


sin  6  =  sin  i  (g^) 

At  aiv  point  along  its  orbit,  the  velocity  conqjonents  of  Mara  are  given 
approximtely  by  the  expressions 

“  Pm  -  2®M  " 

(25) 
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The  above  relations  are  used  to  obtain  a  nev  expression  for  the  charac¬ 
teristic  transfer  velocity#  as  shown  In  Eg,.  (26). 


(26) 

(a  derivation  of  Eg.  (26)  Is  given  in  Appendix  B. ) 

As  Is  to  be  expected#  Eg.  (26)  beecmes  Identical  with  Eg,  (5)  If  we  set 

Sm  -  €  »  &  »  0 

Ihe  steps  necessary  for  the  solution  of  the  present  case  are  not  too 
different  from  the  ones  described  in  the  coplanar  case,  ^e  main  difference 
arises  in  the  computation  of  the  time  of  transfer, 

At  the  time  of  departure#  the  InltleuL  angular  position  of  the  Earth,  ro, 
and  of  Mars#  0^,  are  known  quantities.  If#  as  before#  a  transfer  angle  6  Is 
chosen#  one  can  coiqpute  the  mean  anomalies  and  from  Kgs.  (20)  and  (21). 
Similarly#  with  known#  pjj  can  be  obtained  from  Eg.  (I8). 

With  the  above  Information  we  can  now  find  the  length  of  the  ciiord  Joining 
points  P  and  Q  from  the  relation 

-  1  +  -  2Pu  cos  9  (27) 

2  2 

where 

Py  «*  magnitude  of  Pj^  idien  M  ■  in  Eg.  (I8) 
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liavins  now  obtained  the  valiies  of  t  and  C,  ve  can  continue  with  the  com- 
putatlonal  procedure  outlined  in  the  previous  section.  Eq.mtlons  (llf)  and 
(s4)  ^ow  that  for  given  fixed  values  of  at,  0,  and  1>  after  account  is  ta4en 
of  Eq.  (27)y  the  angles  &  and  e  will  depend  only  on  the  chord  C«  Erom  this 
we  conclude  that  the  condition  of  optimality,  Eq.  (l2).  Is  also  applicable 
in  the  present  easoj  provided  Eq.  (22)  is  employed  to  conpute  the  tins 
appearing  in  Eq.  (9), 

Hie  e:qpreaalon  for  the  inclination  of  the  departure  velocity  vector 
remEd.ns  unaffected  by  the  eccentricity  ej^  and  can  be  used  in  Its  Initial  form. 

It  was  noted  that  for  the  two-dimensional  circular  case  the  only  input 
Infonaation  needed  was  the  initial  longitude  difference  of  the  two  planets. 

In  the  present  case,  their  actual  positions  in  orbit,  w  and  0^)  have  to  be 
specified.  In  all  other  respects  the  method  of  carrying  out  the  solution 
remains  undianged. 

Figure  9  shoirs  a  plot  of  vs  time  for  the  three-dimensional  model  of 
the  planeteoy  system.  In  es<di  case,  departure  was  assumed  to  take  place  on 
October  1^  i960. 

It  can  be  seen  that  the  curves  consist  of  two  separate,  parabollcally 
shaped  branches.  One  branch  corresponds  to  values  of  the  travel  angle  6  < 

180  deg^  while  the  other  refers  to  9  >  160  deg.  Due  to  the  inherent  geooietzy 
of  inclined  planes,  txanafere  throu^  angles  in  the  immediate  vicinity  of 
180  deg  are«  in  general*  laposelble  except  for  tranefer  trajeetorlee  that 
orlglmte  at  the  nodal  line. 

I’ve  reglone  of  minloal  transfer  velocity  are  now  found  for  each  seleoted 
daparture  date,  separated  from  each  other  by  the  narrow  belt  6  180  deg. 

Ihe  value  of  on  which  the  optimization  prooedure  outlined  before 


Fig.9~Three-dimensional  characteristic -velocity  variation 
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vlU  converge,  will  depend  on  the  initial  value  of  0  selected  to  start  the 
iteration.  The  for  long-period  transfers  will  he  foiaid  if  one  starts 

■OtM  iteration  with  a  value  of  0  >  180  deg,  vdille  the  short  transfer  time 
Vjjjj  is  reached  upon  starting  out  with  a  9  <  l8o  deg.  An  idea  of  the 
change  in  orientation  of  the  departure  velocity  vector  as  the  time  of  trans¬ 
fer,  as  well  as  the  transfer  angle,  is  increased  can  he  gained  from  Fig.  10. 
Ihis  figure  presents  the  two  angles  S  and  plotted  against  transfer  time  t. 


Time  of  travel,  t  (days) 


_L 


-L 


J_ 


J 


100  l‘»0  180  220  260 

Transfer  ongle,  9  (deg) 


300  340 


Fig. 10 — Variation  of  8  and  y.  with  flight  time 
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Appendix  A 
NUMERICAL  EXAMPLE 

A  tnsaerleal  exaaple  Is  vorked  out  In  detail.  In  this  appendix  to  Illus¬ 
trate  the  general  method  of  solution. 

For  the  sake  of  convenience^  a  two-dinensiorval  circular  planetsixy  model 
has  been  selected,  with  the  following  input  ilp.ta: 

9  n  Initial  longitude  difference  »  30  deg 
(Mars  leading  the  Earth) 

$  >  transfer  angle  »  IkO  deg 

■  1.0  «  radius  at  point  of  departure  P 

Pg  ■  1.523  "  radius  at  point  of  arrival  Q 

From  the  above  data,  the  chord  joining  points  P  and  Q  is  coirputed  to  be 

C  -  1  +  1.523®  -  2  •  1.523  cos  140°  ■  2.3776 


Die  exanple  chosen  here  happens  to  have  a  vsdue  of  the  transfer  angle 
9  <  ir,  so  that,  in  accordance  with  the  discussion  presented  in  connection 
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with  Fig.  3,  the  appropriate  es^iresaion  for  the  time  of  transit  would  have 
to  he  selected  from  Kqs.  (la)  or  (ih)  (t^^  and  selection  can 
he  narrowed  down  further  if  we  detenalne  whether  Is  greater  or  smaller 
than  iaie  value  of  t  at  A  a  at  'Aloh  point  «  "^EF^Q* 


‘PS 


^^3/2  q)  _  (P  _  ain  $) 


For  Idle  present  case 


1  +  p  +  C 

S  « - i -  a  2.4503 


S  -  C  a  0.0727 


thus 


P_  ■  2  sin 
Xfi 


sin  “  0 


sin  Pjjj  -  0.3393 


T  -  (1.22515) 

SI 


3/2 


3.1416  -  0.34-62  +  0.3393  -  4.251 


Now  that  we  have  determined  that  we  conclude,  with  the  aid  of 

Fig.  3,  that  the  appropriate  time  expression  to  bo  used  in  the  present  case 
is  the  tpQ  relation;  the  vacant  focus  of  the  proper  transfer  ellipse  will  be 
located  on  the  branch  of  the  ciarve  of  Fig.  2. 

What  remains  to  he  done  now  la  to  solve  from  the  expression  for  Tp^, 

E(i.  (la),  for  the  correct  value  of  A  which  satisfies  the  condition 
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This  vill  be  aecooplished  ty  a  tiial-and-error  procedure  •vdiloh  can  be 
suitably  adapted  to  locate  the  proper  root  very  quickly. 

Let  us  start  out  vith  a  first  guess  for  A. 

Queas  1;  A  *»  1. 32 
a3/2  «  1.5167 

a  »  2  . 2  sin"^  O.963I+  »  2.5988 

p  «  2  sln"^^-^^»  2  sin"^  O.I6595  o  O.3334 
sin  a  “  0,51654  sin  p  •  0,32726 


and 


T 


1.5167 


2.5988 


0.51654 


0.3334  +  0.32726 


3.1489 


Since  this  value  is  smaller  than  the  desired  am,  ve  vould  have  to  take 

a  smaller  valw  of  A  as  a  second  guess.  In  order  to  cone  up  with  an  intelll> 

gent  second  guess  for  A,  ve  can  make  use  of  the  shape  of  the  t  vs  A  curve 

of  Fig.  3,  in  particular  that  portion  in  the  neighborhood  of  the  point  (t.  A  ). 

nr  n 

If  ve  drev  a  curve  between  the  latter  point  and  the  point  {3. 1489,  I.32)  as 
shown  in  Fig.  11,  such  that  dr/dA  m  „  at  ve  find  that  at  a  value  of 
T  ■  3.6,  A  a-  1.25. 

Accordingly  ve  let  our  second  guess  be 
Guess  2;  A  ■  1.25 
a3/2  «  1.3975 

a  -  2  sln"^  0.99001  -  2.8588 

p  -  2  -  sln"^  O.I7O53  -  0. 34276 


sin  a  ■  0.2867 


sin  p  ■  0.33608 


Transfer  time 
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Semimajor  axis,  A 


Fig.ll— Approximated  t  vs  A  curve  for  graphical  iteration 


-39" 


and 


T 


1.3975 


2.8588 


0.2867 


0.34276  +  0. 


3.5851 


Plotting  a  corrected  curve  In  Fig.  11  throu^  tke  three  known  points, 
ve  find  that  a  third  guess  should  be  around  A  >  1.2485. 

Quess  3;  A  ■  1.9485 

Bepeating  the  calcxilations  perfoxned  before,  ve  arrive  at  a  value  for 
time  T  of  3.613.  To  the  degree  of  accursuqr  needed,  a  satisfactory  solution 
can  be  taken  to  be  a  value  of  A  a  1.2487. 

QSie  total  characteristic  transfer  velocity  can  now  be  counted  from 
Bq.  (5).  First  ve  shall  ccs^pute  the  constants  and  Kj^  appearing  there: 


K 


2Uj> 

3"^ 


where  v„  and  v„  _ _  are  the  surface  escape  velocities  for  Earth  and  Mars, 

esc  n  esc  ' 

respectively. 


OJie  magnitude  of  the  senilatus  rectum  L,  calculated  from  Eq.  (2a),  is 
found  to  be  L  ■  I.I87. 

With  the  above  values,  Eq.  (5)  gives 


''or,  ■  {3  -  ISW  ■  •* 


1/2 


1:2k  3  -  -  ^ySif  *  1-523  •  0.0861}' 


1/2 


0.4020  +  0.1944  -  0.5964 


-1*0. 


Slziee  'f;he  two  neaibers  of  tbs  above  sum  supply  the  velocity  requlxeients 
at  depeurture  and  arrival,  re  spsetl valy,  we  find  in  dltoensional  form 

-  0.59611-  •  97#  680  -  58,256  ft/sec 
V„  -  0.40S0  •  97>68o  «  39,267  ft/sec 

CHi 

V_  «  0.19l^l^  •  97,680  -  18,989  ft/sec 

CHg 

Ihe  angle  7^  \diloh  the  departure  velocity  vector  subtends  with  the  Baith's 
oztltal  path  la  found  to  be 

7^  -  tavT^^ -  m  tan'^  O.IOI23  -  sV’ 

Ihls  completes  the  Illustrative  computation  for  the  two-dimensional 
exasQila  chosen.  Ihe  same  basic  scheme  can,  of  course,  be  followed  also  in 
the  more  eonplioated  ease  of  three-dimensional  clrcular-to-elllptie  transfers, 
provided  the  tins  Is  properly  evaluated  and  the  radius  Pu  has  been  de- 

O  Mg 

termlned  In  accordance  with  the  method  outlined  previously. 
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Appendix  B 

CRARACTERianC  TRANSFER  VELOCm  IM  IHBKB  DIMENSIONS 


ilhe  velocity  coogionents  of  a  'bo<^  moving  along  a  conic  aection,  pjrojeeted 
onto  the  radial  and  transverse  directions,  were  given  hy  Eqs.  (2.142)  and  (2.1.3) 
of  Ref.  1.  Applying  these  e:Q)resslons  to  the  planet  Mars  ve  have 


(B.1) 


In  the  above  relations,  denotes  the  Instantaneous  radial  distance  of 
>kr.  *11.  Tj  tort.,  th.  maiu.  ot  the  ortlt  *lct,  to  flr.t-or4.r 

terms  in  e^,  eq;uals  the  seoimajor  axis  a^^.  !nie  tens  Ij^  Is  the  semllatus 
rectum  of  the  Martian  orbit. 

If  one  enqilayB  the  relations 


1  >  cos 


M 


(B.2) 


differentiates  the  flrat  to  obtain  and  substitutes  both  Into  the  v^ 

expression  of  £q.  (B.1),  one  finds  that  for  snail  values  of  Ou 

H 


sin  M 


•* 

'75  (l..„.o.  M) 


(B.3) 
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Denoting  ty  the  velocity  of  the  vAicle  relative  to  a  nonattracting 
planet  located  at  point  Q,  ve  can  vrlte 


v2 


■  ’''J 


*(M); 


H 


8ln®6 


(re)jj  cos  e  -  v. 


m 


(B.4) 


\diere  (r)jj  and  (rd)y  "the  heliocentric  valociiy  oon5ponente  of  th«  vehicle 
In  its  own  plane  of  notion,  evaluated  at  the  point  of  intersection  of  the 
transfer  plane  with  the  MEortlan  orbit,  and  are  given  by 


(B.5) 


The  symbols  1  and  a  refer  respectively  to  the  senllatus  rectum  and  semi- 
major  axis  of  the  transfer  ellipse*  ISquatlons  (B<3)  and  (B.3)  can  now  be 
substituted  into  Eq.  (B.4).  After  e>;pandlng  some  of  the  terns,  retaining 
only  linear  terms  in  ej^  and  combining  them  appropriately,  one  is  led  eventu¬ 
ally  to  the  following  nondlmensionalized  relation  for  vf- 


is  ■  fc  ^3  -  ?  -  ( 1  -Jh  '»•  ‘ ) 


(B.6) 


-  2e„  sin 


f 


Pg  A 


For  the  departure  portion  of  the  transfer,  the  angle  6  tabes  the  place 
of  the  angle  e  used  above.  The  relative  departure  velocl'^  Increment  is 

■  A  ■  ^  cos  b|  (B.7) 
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the  relation  ataown  In  Sq.  (26) 
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